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Abstract 

A consistent approach (based on QFT) to neutral fermion creation (due to their 
[ magnetic moments) in strong inhomogeneous magnetic fields is considered. In particu- 

lar, it is demonstrated that in specific cases, the problem can be technically reduced to 
i-^ the problem of charge particle creation by an electric step. We demonstrate that the 

latter problem can be solved consistently in the framework of QFT. All the details of 
such a construction will be described by us in a next publication. Here, however, using 
^ ' some particular results of the general construction and the above mentioned analogy, 

\Q , we calculate neutral particle creation from the vacuum by a linearly growing magnetic 

C-» ■ field. 

1 Introduction 

Usually, particle creation from the vacuum by strong electromagnetic fields is associated 
with creation of charged particles by strong electric-like fields. Acting on virtual charged 
particles an electric-like field can produce a work and materialize them on the mass shell 
as real particles. Nevertheless, if a neutral particle has an anomalous magnetic moment, a 
inhomogeneous magnetic field acting on such a particle, can also change its kinetic energy 
(produce a work). This mechanism can provide neutral particle creation from the vacuum 
by strong inhomogeneous magnetic fields. In this respect, one can speak about two can- 
didates (among the known elementary particles), a neutron and a neutrino. It is known 
that the neutron has negative magnetic moment given by fi n = —1.9130427(5)/%, where 
fi N is the nuclear magneton, [i N = e/2rnjv- It is also possible that neutrinos have magnetic 
moments, for a review see [Tjr2].The recent experimental constraints on the neutrino mag- 
netic moments are in the range fi v < 3.2 x 10 -11 /i s 3 , where [i B — e/2m e is the Bohr 
magneton. Astrophysical constraints on the magnetic moment of the Dirac neutrino can be 
even stronger [3]. Note that in order to satisfy m v < leV, the theory argues that a more 
natural scale for the Dirac neutrino would be \i v < 10~ 14 /x B [5]. The effect under discussion 
can be observed in inhomogeneous magnetic fields that have to be very strong in a certain 
domain. Such fields can exist in the nature. It has been suggested that magnetic fields of 
order 10 15 - 10 16 G or stronger, up to 10 18 G, can be probably generated during a supernova 
explosion or in a vicinity of the special group of neutron stars know as magnetars, see, for 
example [B] . For magnetar cores made of quark matter the interior field can be estimated to 
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reach the values B ~ 10 20 G [7J. The possibility to create a strong quasi- uniform magnetic 
field with the strength of the hadronic scale B ~ 10 19 G, or even higher in heavy-ion colli- 
sions at RHIC and LHC, when the matter in the central region is presumably in the QGP 
phase, is recently shown [5]. Superconducting cosmic string - if they exist - could generate 
fields more then 10 30 G in their vicinities [5]. 

Recently, the Schwinger effective action approach [TO) was formally applied to calculate 
the probability for the vacuum to remain a vacuum in a linearly growing magnetic field for 
neutral fermions of spin 1/2 with anomalous magnetic moment. The same problem in 2 + 1 
dimensions was considered in [TT] , in 3 + 1 dimensions in [T2] . It is difficult to accept results 
presented in [12j . which, in particular, admits neutral particle creation in a homogeneous 
magnetic field. This means that formal calculations a la Schwinger, without any theoretical 
justification based on quantum field theory (QFT), can lead to mistakes. Results of the 
work [11] seem to be reasonable, but essentially use specific of gamma matrices in 2 + 1 
dimensions, and cannot provide complete description of the effect. 

It should be noted that until now a consistent description (based on QFT) of particle 
creation (due to their magnetic moments) in strong inhomogeneous magnetic fields was 
unknown. To provide such a description is a part of the present article. In the framework 
of the developed approach, we demonstrate that in a specific cases, the problem can be 
technically reduced to the problem of charge particle creation by an electric field given 
by a step scalar potential. We analyze once again the latter problem, using general QFT 
formulation of particle creation by potential steps presented recently, and formulating for 
the first time its solution from the view point of QFT. As a final result, we calculate neutral 
particle creation from the vacuum by a linearly growing magnetic field. 

2 Dirac-Pauli equation with a constant magnetic field 

In 3 + 1 dimensions (dim.) and in the relativistic quantum mechanics neutral fermions of 
spin 1/2 and mass m with anomalous magnetic moment fi (without electric dipole moment) 
in an external electromagnetic field F\ v are described the Dirac-Pauli equation, see [TH [15] . 
Such an equation has the forrrQ: 

(Vpa - m - ^a Xv F x ^j $ (x) = , 

p v = id v , <7 A " = |[7 A ,7l, (1) 

where F\ u (x) is the field tensor, ip (x) is a four spinor, x = (x° =t, r), r = (x,y, z), and 

l" = (7°j7) are Dirac matrices. 

Let the external field is a constant nonuniform magnetic field B that is directed along 
the z-axis and depends on the coordinate y only, B (y) = (0,0, B z (y)) such that the only 
nonzero components of the stress tensor are F21 (y) = — -F12 (y) = B z (y). In addition, we 
suppose that B z (y) takes constant values as y — > ±00, such that 

d y B z (y) 0. 

Moreover, we suppose that for y < yh (the range Sl = (— 00, y^]) and y > yn (the range 
Sr. = [yn, 00)) the field B z (y) is already uniform and its values are B z (y) = B z (—00) and 
B z (y) = B z (+00), respectively. Thus, the magnetic field under consideration is constant 
uniform (or zero) at spatial infinities and, in fact, represents either a potential barrier or 

1 Here we are using the natural system of units h = c = 1. 
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step for the magnetic moment fi. With such an external field, equation ([1} takes the form: 
id iP (t, r) = Hr/> (i, r) , H = ( 7 ° 7 3 P 3 + 7 °£ 2 II 2 



it 



Sz7P-L + m£ 2 - fiB z (y) , p ± = (p 1 ,p 2 ,0) 



(2) 



In the case under consideration, the operators p , p , p 3 , and IT Z are mutually commuting 
integrals of motion (all these operators commute with the Hamiltonian H). The integral 
of motion Tl z is generalization of z-component of a spin polarization tensor for a uniform 
magnetic field, see [T5] . 

It is useful to use an additional spin operator R, which is also an integral of motion 
commuting with the previous ones, 



r = hu: 



1 



ffl- 1 



-1/2 



(3) 



A complete set of solutions of equation @ can be written in the form 

i>n (*! r ) = ex P {-Wot + ip x x + ip z z) ip n (y) , (4) 
where ip n (y) are eigenvectors of the equations 

suj\J\ + {p z /uj) 2 Rip n (y) = p ip n (y) => Rip n (y) = sip n (y) , p = + (Pz/uf, 



R 



-1/2 



1 + { Pz /uj) 2 \ "'" ( S7 7 3 p z /^ + 7 S z ) 
n 2 (p x , y) - suj ip n (y) = 0, s = ±1, 

n* (p x ,y) = - m-Sz (y) > ^ z = s 2 ( 7 Vx + i 2 P 2 ) + m^ z , 



(5) 



(6) 



and n — (p x ,p z ,u,s) is the set of quantum numbers from a complete set of numbers that 
will be specified in which follows. Choosing ip n (y) as 

iPn {y) = \{l + sR)$(y), 

where $ (y) is an arbitrary spinor, we obey eq. ([5]). The real continuous quantum number uj 
can be positive and negative and determines the transversal part of full energy, uj 2 = pg—p 2 , 
that is, it determines the full energy of particle moving on xy plane. 
Then solutions of eq. © can be represented as 



1> n (v) = g (! + sR ) + » B * (y) + «w] 0„ (y) , 
where the spinors <p n (y) satisfy the following equation: 

{-d 2 + m 2 +p 2 x - pr?d v B 9 (y) - [uj + s»B z (y)] 2 } 4> n (y) = 0. 
It is convenient to represent the spinor <p n (y) in the form 

0™ (y) = v n ,x {v) + w' 1 ) w ' 



(7) 



(8) 



(9) 



where it is selected that either x = +l° r X = ~ 1, u is an arbitrary constant spinor, and 
scalar functions ip n x (y) are solutions of the equation 



{~ d y + m2 +Pl + iX^dyBz (y) - [w + s\lB z (y)] 2 } ip n x (y) 



= 0. 



(10) 
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In what follows, we suppose that v is normalized as v^v — 1. In addition, vv^ is the identity 
4x4 matrix, vv^ = I. Thus, the spinor structure of solutions ([7]) is defined completely. One 
can easily verify that solutions ([7]) that differ by values of \ only, are linearly dependent, this 
is an effect which the projection operator [. . .] in the representation (J7J produces. Because 
of this, it is enough to work with solutions corresponding to one of possible two X- That 
is why, the upperscript x will sometimes disappear from solutions, but in such cases it is 
supposed that x 1S fixed in a certain way, the same for all solutions under consideration. 

Using the freedom inherent in the solutions of eq. (|10[) . we construct two (in general dif- 
ferent) sets {<;if> n (^ r )} an d {^V'n (^) r )} 01 independent solutions, £ — ±, satisfying specific 
boundary conditions as y — > — oo or y — > +oo. The first set contains states c,tp n (t,r) with 
definite real values p L of y-component of the momentum, such that £ defines the sign of the 
momentum, 

- id y C V>„ (t, r) = p L c ip n (t, r) , C = sgn p L , y -> -oo. (11) 

The second set contains states ^ip n (t,r) with definite real values p R of the y-component of 
the momentum, and again £ defines the sign of the momentum, 

- id y (t, r) = p R ty„ {t, r) , C = sgn p R 7 y -> +oo. (12) 

We are interested in the nondecaying solutions of eq. ([TU)) as y —> ±oo. In this case both p L 
and p R are real. We believe that for any given quantum numbers n both set {(ip n (A r )} 
and {''ipn (^ r )} represent complete set of nondecaying solutions. In fact this is the above 
mentioned supposition about the form of the field B z (jj). 

It should be noted that the time independence of the magnetic field under consideration 
is an idealization. In fact, it is supposed that a field inhomogeneity was switched on in a 
time instant U n then acts as the constant field during a large time T, and was switched off 
in a time instant t out = U n + T, and one can ignore effects of its switching on and off. This is 
a kind of the regularization, which could be, under certain conditions, replaced by periodic 
in t boundary conditions. Namely, by analogy with periodic boundary conditions in space, 
that are usually imposed as the volume regularization, here we impose periodic (with the 
period T) boundary condition in time t. Thus, we consider a theory in a big 3-dimensional 
space-time box that has a volume V y — TS XZ , S xz — L x x L Zl where L x , L z , and T are 
macroscopically large, L x , L z — ¥ oo and T — > oo. 

It is convenient to use the inner product on the time-like hyperplane y = const, which 
has the form 

(V,V0 V = / ^(i,r)7 7V ' (t,v)dtdxdz. (13) 

JVy 

The integration in (|T3| is fulfilled in the limits from —L x /2 to +L x /2, —L y /2 to +L y /2, 
and from —Tj2 to +T/2 in the time t. It is supposed that all the functions ip are periodic 
under translations from one box to another one. Under these suppositions, inner product 
(|T3|) does not depend on y. We note that the quantity (fl~3")) for t// = ip represents the particle 
current via the hyperplane y = const. 

By using the inner product (|13|) . we obtain: 

Thus, the current density in y-direction in the state ip n (t, r) is 

in = ^ n (y)i°i 2 t n (y)- (14) 

Using the structure ([7]), we rewrite the combination ijr n (y) 7°7 2 '0n (y) as follows 
< x (V) tr | ( 7 V + fdyl 2 ) - mZ z - [iB z (y) - slj] 7 °7^ [1 + sR] 



(7 X Pa 



i^vi 2 ) - m ^z - v B z (y) - su - (1 + xn 1 ) I f'n, x (y) 
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where tr{...} is the trace in the space of 4 x 4 matrices. Calculating this trace, we obtain 

€ (y) 7V< (V) = (l + (PzM 2 ) < x (V) - ^y) + s»B z (y) + s X it y ) f' n , x (y) ■ 

(15) 

As was already mentioned, we supposed that B z (y) tends to some constants values as 
y — > ±oo. Let us suppose for definiteness sake that the derivative d y B z (y) has a definite 
sign, let say d y B z (y) > 0, Vy, and let fi < 0. Note that there are no bound states in this 
case. To simplify the consideration, we also suppose that 

U = U R -U L > 0, U L = -nB z (-co) < 0, U R = -nB z (+oo) > 0. 

For asymptotic (as \y\ — > oo) states with real values p L and/or p R , we have 

CVn, x (y) = C-^exp (ip L y) , VL x) (*) = C AAcxp (ip fl y) , (16) 

respectively, where fjV and vV are normalization factors. We introduce the notation 



E S (L/R) = n s (L/R)Jl + \p z /TT s (L/R)] 



l2 



7T S (L/i?) = W - SU L/R , 71^ = V^l + m2 

and note the following relations 

7T S (L) = TT S (R) + sU (17) 

and 

(p L f = [E s (L)] 2 - nl-pl {p R f = [E. (R)] 2 - nl-pi (18) 

where the last hold due to eq. (H~U|) . We see that \E a (L)\ and \E S (R)\ are the asymptotic 
values of kinetic energy, while \ir s (L)\ and 1 7r s (R)\ are the asymptotic values of its transversal 
part, respectively. 

Then, using the asymptotic conditions (jTTJ) , (|T2"j). and the result (fT5")) . we can subject the 
introduced sets {(ip n (t,r)} and {^ip n (^ r )} to the following orthonormality conditions 

(cVW V\i') x = C»72A,e"W; ( C ^ nl C ' V>„') = C^iA.CW', (19) 

where 

»7l = sgn7r s (L) , t?^ = sgn7r s (i?) . 

In deriving (IT51) . it was taken into account that for asymptotic (as \y\ — > oo) states with real 
values p L and p R , the relations 

\n s (L)\ > \p L \, \tt s (R)\ > \p R \ 

hold due ([TBI , respectively. That is why the sign of the quantity (fTS"]) with the operator 
7r s (L/R) + s\i d y is due to the sign of the n s (L/R). The normalization factors in (fl"6")) 
are as follows 

C Af = C CY, ttf=tCY, Y=(l + ( Pz /u J ) 2 y /i V- 1 ? 2 , 

C C = [21^11^(1,)- ^0 - 1/2 , ^=[2|^||7r s (i?)- SX p R |]" 1/2 - (20) 

In the limit of infinite volume of the normalization (continuous momenta po, p x , and p z ) one 
has to substitute 5 n , n ' in normalization conditions flT9"|) by 8 StS >S (po — p' ) S (p x — p' x ) 5 (p z — p' z 
In this case, Vy 1/2 -> (2tt)~ 3/2 in eqs. {2D]). 
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For any given quantum numbers n, both set {(tp n (t, r)} and Vn (*j r )} represent com- 
plete set of nondecaying solutions of equation ([2]). Their mutual decompositions have the 
form: 



Vl ^ n (t, r) = (t, r) <? ( + |C ) - _ V„ (t, r) 3 (_ |< ) ; 

VRC^n(t^) = + ^n(*= r )ff( + lc) - r ) .9 (~ k) . 

where the decomposition coefficients g are defined by the relations: 

(A, C 'Vv) = <Ws(c| C ')< .9( C 'k) ^ .9 (c | C ')* ■ 



(21) 



(22) 



Using the orthonormality conditions (1191) . we derive the following relations for the decom- 
position coefficients: 



.9 



.9 



.9 1 



9(c\ + )9( + k)~9(c\-)g(-k) = OiLVRSt,c 

In particular, these relations imply that 



- \ |2 



k(-i + )i=k( + i-)r. | 5 ( + i + )r=M-i-) 



g(+r) g( + l-) 

5-C-I-) s( + l+)' 



(23) 



(24) 



Thus, one can see that all these coefficients can be expressed via only two of them, e.g. via 
g (+ | + ) and g { + \~ ). However, even these coefficients are not completely independent, they 
are related as follows: 

\9(+\-)\ 2 -\9(+\ + )\ 2 = -VlVr- (25) 



3 Creation of neutral fermions 

ft is useful to make a preliminary qualitative analysis of the behavior of particles and an- 
tiparticles in the fields under consideration. It should be noted that here there exist two 
principally different cases, the first one corresponds to U < 2m, whereas the second one 
(we call it creation case, C-case) corresponds to U > 2m. In the first case, there exist only 
a scattering of neutral fermions by the magnetic field without additional particle creation 
from the vacuum. This case can be treated in the framework of one-particle relativistic 
quantum mechanics. The quantum number s gives the spin polarization both for particles 
and antiparticles. Choosing the magnetic moment of the particle as [i, we have the mag- 
netic moment of the antiparticle as — fx. Note that we fix // = — \fi\. Then, according to the 
standard particle-antiparticle identification of wave functions, asymptotic kinetic energy (at 
y ±oo ) of the particle moving on xy plane is ir s (L/R) > 0, while it is — tt s (L/R) > 
for the antiparticle. One can see from eq. (IT71) that the particle potential energy s B z (y) 
decreases along the axis y for s = — 1 and increases for s = +1. At the same time, the an- 
tiparticle potential energy — s \fi\ B z (y) increase along y for s = — 1 and decrease for s = +1. 
That means that the field B z (y) accelerates particles with s — —1 and antiparticle with 
s = +1 along the axis y. Respectively, antiparticles with s = —1 and particles with s = +1 
are accelerated by the field in the opposite direction. The same observation holds in the 
case U > 2m. 

We note that real particles are described by some localized in the space-time wave pack- 
ets, such that we have to study the motion of such packets in the external field (obviously, it 
is enough to speak about a localization in the y-direction). Let us denote by Si n t the region, 
where the magnetic field is inhomogeneous. In the region Sl, situated left from S{ n t and in 
the region Sr right from S- ln t, the magnetic field is homogeneous. For big enough differences 



G 



U between the initial and final potential energies, particles and antiparticles with any initial 
kinetic momenta along the axis y get final kinetic momenta directed always as their the 
acceleration by the magnetic field. That is what we have in the case U > 1\Jp x + pi + m 2 
for all partial waves with given p x and p z of a wave packet. Because particles and their 
antiparticles with a given s have opposite directions of the acceleration, there exists a state 
polarization out of the region Si n t- The final particles with s = +1 and antiparticles with 
s = — 1 are situated in the region S l , and final antiparticles with s = + 1 and particles with 
s = — 1 are situated in the region Sr. 

From the physical point of view, there is a similarity between the two cases, one where 
neutral fermions with anomalous magnetic moment are placed in an inhomogeneous mag- 
netic field B z (y) with d y B z (y) > 0, and another one where charged fermions are placed in a 
constant electric field directed along y and given by a scalar potential A (y) . In both cases 
external fields produce a work which implies an acceleration of the corresponding particles in 
y-direction. From the QFT point of view if such a work is greater than 2m (C-case), a parti- 
cle creation from the vacuum is possible. In fact, this analogy allows in both cases formally 
to use the same techniques of calculation. It turns out that the problem of neutral fermion 
creation in strong inhomogeneous magnetic field can be technically reduced to the problem 
of charge particle creation by an electric potential step. Some heuristic exact calculations of 
the particle creation by potential steps in the framework of the relativistic quantum mechan- 
ics were presented by Nikishov [16l [17] , developed then in p~8j , and were used in numerous 
works in the frame work of semiclassical considerations, see reviews [19] [20] . 

In such a way it seems that we could use the known results to find the mean number 
of neutral particle-antiparticle created. However, a closer consideration shows that the 
particle-antiparticle and causal identification of wave functions ^tp n (t, r) and £ip n (t, r) given 
by Nikishov [T6l |T7] does not coincide with one given by Hansen and Ravndal 18 for the C- 
case, see discussion in [21]. Within the WKB approximation this difficulty can be bypassed, 
but the question remains. Trying to resolve this contradiction, we have realized that at that 
time no any justification for quantum mechanical calculations from the QFT point of view 
were elaborated. Such a justification is given in our recent work |22) as a part of general 
QFT formulation of particle creation by potential steps. We follow the principle ideas of 
this work in solving the problem under consideration, taking into account some necessary 
identification. 

In the C-case, there exists a range 2 yjn 2 , + p z < U of momentum p x and p z of the 
fermions, such that particle creation is possible. This case is described by wave functions 
([7]) with quantum numbers of domain f2, where uj,p x , andp z are restricted by the inequalities 



sir s (L) > ir x , stt s (R) < -ir x , 2v^2 +p 2 z < U. (26) 

If to treat this case using identification of wave function by an analogy with one-particle 
scattering theory, there appears an analog of Klein paradox for charged relativistic particles 
in an electric field [23] . This is an indication that one has to use appropriate many-particle 
description given by QFT to treat the problem correctly. 

On the first stage of the canonical quantization of the field ip (t, r) one establishes that 
the corresponding quantum field is the Heisenberg field operator <J/ (i, r) that satisfies the 
equal-time anticommutation relations: 

[* (t, r) , * (t, r')] + = [* (f, r) f , * (t, r') 1 ] = 0, [* (t, r) , V (t, r') f ] = 5 (r - r') (27) 

and the the Dirac-Pauli equation ([2]). The formal expressions for the Hamiltonian T~L of 
the quantized fermion field and the corresponding magnetic momentum operator M. can be 
easily constructed, 

H = J $(t,r) f HV(t,r)dr, M = | / * it, r) f , * (t, r) dr. (28) 
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To perform quantization in terms of particles and antiparticles, we define the following 
inner product 

(tM') t = / 1>Ht,T)il>'(t,T)de (29) 



between two solutions of the the Dirac-Pauli equation on i-const hyperplane. This inner 
product does not depends on the choice of such a hyperplane if the spinors ij) (t, r) obey 
certain boundary conditions that allow one to integrate by parts in (|29[) neglecting boundary 
terms. Since physical states are wave packets that vanish on the remote boundaries, the 
above supposition holds true and the inner product (|29[) is time independent for such states. 
Considering plane waves instead of natural wave packets, one has to impose corresponding 
periodic boundary conditions on the corresponding wave functions and external field to 
maintain the inner product (|29[) time independent. However, in the case under consideration 
the external field with different asymptotics at y — ¥ ±00 cannot be adopted to any periodic 
boundary conditions in y-direction without changing its physical meaning. Here to provide 
time independence of the inner product, one has to modify the inner product itself. This 
modification is applied to the integration over y in the expression (|29p and is described 
below. 

Let ip n (t, r) and i/j' n , (t, r) be wave functions ([7]) and the integral over the variable y in 
the infinite limits be regularized by large positive numbers L\ and L%. Integrating over the 
variables x, y, and using representation @, we obtain 

/L2 
Qdy, 

Q = (<Pn, x (y)Y 1t I + {uj + snB z (y) + sxidyf <p' na (y), (30) 

where the orthogonality for n 7^ n follows as L\, L2 — > 00. 
We represent the regularized integral TZ as follows 

/Vl i-Vr rL 2 

Qdy+ / Qdy+ Qdy, (31) 
Li JyL Jvr 

where only the second term, integral over the region S- m t, depends on the derivative d y B z (y). 
The smoothness of the d y B z (y) allows us to believe that this integral is finite as L\, L2 —> 00. 
The first and the third terms are calculated as integrals over the regions where d y B z (y) = 0. 
Then their values are determined by asymptotics (I16|) in the following form 

/Vl rL 2 
Q L dy, Ur= Qndy, 
Li Jy R 

Ql/r = (^ x (y)y[Tr 2 x + (Tr s (L/R) + SX id y f}cp' ntX (y). (32) 

Ql and Qr are constant then TZl ~ £1 and TZr ~ L2. We see that only TZl and TZr are 
essential in (1311) as L\, L2 — > 00, 

TZ — > K L +K R . 

There exist two independent solutions with given quantum number n of domain tt. In 
spite of the fact that these solutions are obtained in the constant external field we believe 
that they represent asymptotic forms of some unknown solutions of the Dirac-Pauli equation 
with an external field d y B z {t,y) that is switched on and off at t — > ±00 and effects of the 
switching on and off are negligible. Since the inner product (|29]) does not depend on t 
for such solutions, we believe that orthogonal pairs of solutions that describe alternative 
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particle/antiparticle states at the initial and the final time instants remain orthogonal at 
arbitrary time instant. Therefore we have to find out which solutions among those we have 
introduced before are such orthogonal pairs. Taking into account the relations (|2"T|) . one can 
show that 

(cVv, C O t = > ne "< ( 33 ) 
if we assume that L\ and L2 satisfy the relation 



Li 



(R) 



0(1). 



(34) 



Condition (13"4"|) guarantees that the sets {^n (^ r )} an d {^V'n (^ r )} f° r n £ ft correspond 
to alternative physical states. 

Consider the quantities R-l/r (|32p defined by the functions Qip n (x) and ^ip n (x) with 
quantum numbers of domain f2. In this case we attribute the corresponding index £ to these 
quantities as follows: R-l/r —> qRl/r. or R-l/r — > ^Rl/r- Using eqs. (fl6|) and (|20|) and 
retaining only leading in the limit L\, L% 00 terms, we obtain 



:Rl 



Y 2 Li 



P L 



C Ri 



Y 2 U_ 



(R) 



P 1 



(35) 



To calculate the quantities qRr and ^Rl, we use relations (|21[) . Again retaining only leading 
in the limit L\,L 2 — > 00 terms (neglecting in particular oscillating terms) and taking into 
account eqs. (fTB"]) and ([20]) . we find 



(Rr 



Y 2 L 2 



Y 2 L X 



tt s ( R) 



p R 
7T S (L) 



9(C 



v 1 - 



+ He 



a - 



(36) 



Note that ^R L > and ^Rr > (Rl due to \g(+ \ )f > 1. Taking unitarity relations 



(37) 



5]) and condition (|34l) into account, we obtain the following orthonormality relations 

7T S (i?) 



p 1 



One can see that the following symmetry take place: Particles with opposite values of s 
have opposite accelerations, the same is valid for antiparticles. That is why the cases s = +1 
and s = — 1 differs only by opposite directions of all the motions, and respectively by the 
opposite dispositions of all the asymptotic ranges. Probabilities of all the processes are equal 
in both the cases. That is why it is enough to consider only one case, let say s = +1. 

It is supposed that we know the complete set of the solutions of the Dirac-Pauli equa- 
tion, parameterized by a set of quantum numbers n, on the hyperplane t = const. Then 
we can decompose the quantum Heisenberg field operator (t, r) and its Hermitian con- 
jugate ^l* (t,r) in this complete set using the inner product (|29|) . Assuming that both 
set {+ip n (t, r) , + ip n it, r)} and {-ip n (t, r) ,~ ip n it, r)} represent the complete set of non- 
decaying solutions for the domain J7, we introduce the notation \I/ n (i,r) for the compo- 
nent of the quantum field operator that can be expanded via either +ip n (t, r) , + ip n (t, r) 
or -t/j n (t, r) ,~~ ip n (t,r). Operator coefficients in such decompositions do not depend on 
space-time coordinates because both quantum field operators and classical solutions obey 
the same Pauli-Dirac equation. For example, for s = +1, we can decompose the ^ n (t,r) 
and (t, r) as follows 



(t, r) = M- 1/2 [a n (out) (t, r) + ftt ( ou t) (t, r)] , 
tft (t, r) = M- 1/2 k (out) +?pl (t, r) + b n (out) +^ (t, r) 



(38) 
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and 



* n (t, r) =M~ 1/2 [a n (in) (t, r) + ftt (in) _ (i, r)] 
*t (t, r) = AT 1 / 2 Lt (in) -^t (t) r) + 6n (in) (tj r) 



(39) 



In what follows, we interpret all a and 6 as annihilation and all a' and 6^ as creation 
operators; all a and cr as describing particles and b and 6^ as describing antiparticles; all 
the operators labeled by the argument in are in-operators, whereas all the operators labeled 
by the argument out are out-operators. It can be shown that these creation and annihilation 
operators obeying canonical anticommutation relations, 

[a n (in),<4(in)] + = [a„(out), aj.(out)]+ = [6„(in), &j.(in)] + = [6 n (out),&£(out)] + = S nk , 
[a„(m),6 n (m)] + = [a„(in), &^(in)] + = [a„(out), 6„(out)] + = [a n (out), 6+(out)]+ = 0(40) 

due to relation (|27[) . In such an interpretation, the in-vacuum |0, in) and out-vacuum |0, out) 
are defined by conditions, 



a n (in) |0, in) = b n (in) |0, in) = 0, Vn; 

a n (out) |0, out) = b n (out) |0, out) = 0, Vn. 



Let consider magnetic momentum operator, 



y n {t,r)\^ n (t,r) 



dr 



(41) 



(42) 



and operator of kinetic energy of the quantum fermion field 
*o (i,r) f 



njkin 

Tin — 



n 2 



»B Z (y) 



1 



p-' 



U z + fiB z (y) 



*o (t,r) dr-ff° 



(43) 



1/ fein 



0,in) is the constant 



in domain tt, where (t, r) = Y^nen *™ r ) an< ^ = \0, in 
term corresponding to the energy of vacuum fluctuations. Using relations (|55|) - (|3"7|) . (|2ip. 
and (|25p. one can represent these operators in equivalent diagonal forms as follows 



Ms 



n 



kin 



t 1 K( m )an(in) - bl(m)b n (m)] 
nefl 

\i ^2 [a* (out)a„(out) - 6* (out) &„ (out)] ; 

[ ~£ n al (in) a n (in) - (in)& n (m)] 
^ [ + £„aj, (out)a„(out) - + £„&+(out)b n (out)] 



(44) 



where 



c £„=£ +1 (i?) + -(7|g(+|-)| 



•£„ - £ +1 (L) - -C/ | 5 (+ 



see details in [35] . We suppose that in the external field under consideration we have 

c £„>0, c £„<0, (45) 

so that the signs of energies ^£ n and c,£ n are determined by the signs of 7r + i {R/L). In 
known solvable cases inequalities (jUJ) hold true, for example, see [HI [13 [21]. Thus, the 
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operator Hq" is positive defined. This fact provides a consistent quantization in terms of 
particles and antiparticles in the domain Q. 

Kinetic energy must be positive for any wave packets of both particle and antiparticle. 
That is why particle wave packets are situated in the region Sl and antiparticle wave packets 
are situated in the region Sr, that is, there is the total reflection from S[ n t both for particles 
and antiparticles. This is consistent with the physical meaning. Note that the expressions 

(<^», C 'Vv) x and (-1)( C ^„ 

'f'^n'Jx' S iven by Eq- (H3), are the probability currents of 
particles and antiparticles through the surface y = const, respectively. The particle and 
antiparticle currents are positive for £ = — 1 and negative for ( — +1. Thus, we see that for 
s = +1 the functions + tp n (t,r) and + ip n (t,r) describe outgoing particles and antiparticles, 
while the functions ~ip n (t,r) and -ijj n (t,r) describe incoming particles and antiparticles, 
respectively. The particle- antiparticle and causal identification of wave functions (J7|) is 
unique in framework of QFT. 

The vacuum corresponds to the absence of incoming particles and antiparticles. In such 
a case the presence of outgoing particles and antiparticles indicates particle creation from 
the vacuum. The effect of particle creation implies constant currents of outgoing particles 
and antiparticles. These currents are equal in the regions 5*^ and Sr. 

Then taking into account eqs. (|38[) and (|39l) . we obtain direct and inverse linear canonical 
transformations between the in and out creation and annihilation operators (Bogolyubov 
transformation) : 

a n (out) = g ( - | + ) 1 g( + \+) a n (in) - g (_ | + ) 6* (in) , 

&t(out) = g (" l+T 1 a n (m)+g^\ + y 1 g( + \+)bi(m); 

a n (in) = g (+ |_ ) _1 g (" |_ ) a n (out) + g (+ |" ) _1 6^ (out). 

bl (in) = -g(+\_y 1 a n (out) + g (+ |- y 1 g (_ |" ) ftt ( out ), (46) 

These transformations are similar to that used by Nikishov in the problem of charge particle 
scattering on electric step [THl [T7] . 

By the help of transformations (|46p. we calculate differential mean number of created 
particles and antiparticles 



N r y = (O,in|at(out)a n (out)|0,in) = \g(-\ + )\ , (47) 
Nt } = (0, in| bl (out) b n (out) |0, in) = | 5 (+ | - ) p 2 . 
Relations imply the equality 

Af(+) = M(-) = N 

which allows us to treat N n as differential mean number of created pairs. The total number 
N of created pairs is the sum 

N=Y,N n - (48) 

The elementary relative probability amplitudes of particle creation, annihilation and 
scattering are defined as follows 

w (+\+) n 'n = c^ 1 (0,out \an> (out) a], (in) I 0, in), 
w (-h)nn' = c v out \ b n> (out) bl (in) | 0, in) , 
w (°l " +)„„/ = c-^out |fet (in)4,(in)| 0,in) , 

w (+ - |Q) /fl = c-^O, out \a n , (out) b n (out)| 0, in) , (49) 
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where c v is the vacuum-to-vacuum transition amplitude, c v — (0, out|0, in) . Amplitudes 
(|49|) are diagonal 

w i+\+)n>n = S n,n'Wn (+1 + ) , W (-]-)„„/ = S n<n 'W n (-|-) , 

w (0| - +) nn , = S n , n ,w n (0| -+),«,(+- |0)„,„ = 5 n , n ,w n (+ - |0) (50) 
and can be expressed via the coefficients g |f J as follows: 

w n (+\+) = g{ + \-)g{-\-r 1 =g{ + \-)g{+\ + r\ 
™ n = g(-\+)g(-\-y 1 =g(-\ + )g(+\ + y\ 
w n (+-\0) = g( + \ + y\ w n (0\-+) = -g{-\-y 1 , (51) 



where transformations (|46[) are used. 

One can express probabilities of particle scattering and a pair creation when n £ f2 , and 
the probability for the vacuum to remain a vacuum via the differential mean numbers N n 
as follows 

P(+\+)n,n> = | < 0,out\a n (out)al,(in)\0,in > | 2 = 5 n>n > -—- P v , 
P(— + \0) n ,n' = I < 0,out\h n (out)a n '(out)\0, in > | 2 = 8 n , n > 1 _ n M Pv , 



P 



r | C-y | 



expJ ^ln(l-iV„)l , (52) 



see details in |22) . The probabilities for an antiparticle scattering and a pair annihilation 
are described by the same expressions P(+\+) and P( — h |0), respectively. 



4 Quasilinear magnetic field 

Here, we consider a specific case of inhomogeneous magnetic field, namely a field linearly 
growing on an interval L y . More exactly, the field of the form 

f B , y < 
B z (y) = I B + B'y, y G [0, L y ] , 
[ B + B'L y , y > L y 

where B' > and Bo = — B'L y /2. Let us call such field quasilinear magnetic field. Consider 
the case given by the condition 

y/\liB'\Ly » max {l, m/VI/^'l} , (53) 

which implies that there is a particle creation in a wide enough range f2 of momenta given 
by condition (|26p . One can demonstrate, similar to the case considered in [3S], see also 
[26] , that leading contributions to the differential mean numbers N n of created pairs do not 
depend on L y in the limit L y — > oo. That is why, it is enough to consider the case of linearly 
growing magnetic field. Equation (fTOj) in the latter field for the function ip n x (y) given by 
(GO) can be written as 



t = v/M^ 7 ^ + ( | M | B'y 1 (|Ml Bo w)J , A = (54) 
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Solutions of this equation, obeying boundary condition (fTTj) and (fT2|) . have the forms 



l<p n>x (y) = M X D- U - X [±(1 + , ±<P n , x (V) = ^*A, [±(1 - , (55) 

where D v (z) are Weber parabolic cylinder (WPC) functions, v = — (iX + 1 + x) /2. By 
the help of an asymptotic expansion of WPC-functions, one can verify the validity of the 
boundary condition (fTTj) and (TT2"1) . Using solutions (1551) . we construct the sets {^ip n (t,r)} 
and (*) r )} °f solutions of the Dirac-Pauli equation. 

The obtained form of solutions formally coincide with the one found in [THl [T7J [25] for 
the case of charged particle creation by constant uniform electric field, compare with |25j . 
Note that our identification of wave functions is in agreement with one given by Nikishov 
for such a special case. This allows us to use these calculations to find differential mean 
numbers of created pairs given by eq. (|47[) . 

As a result, we obtain 

N n = e~ nX (56) 
in the limit y/\p,B'\L y ^> K if the iq and p z satisfies the following condition 



M < w max , \p z \ < w max , w max = \fiB'\ L y /2 - y/\fiB'\K, 

where K is a given arbitrary number K ^> max |l,m/ ^/|/i-B'||. Following the idea of finite 

work regularization presented in |25] , one can show that an exact expression for N n is rapidly 
decreasing as \u)\ — > oo due to the finite work of this field, \uB'\L y , that is, cj max is an 
effective maximum value of quantum number |o>| for the quasilinear field under consideration. 
The maximum value for \p z \ of range f2 follows from condition (|26[) . One can check that 
the mean numbers do not depend on the sign of uB' and on the spin polarization s. Note, 
however, that unlike the case of particle creation due to the electric potential step, the neutral 
particle (antiparticle) created with different s form fluxes aimed in opposite directions. The 
leading approximation given by expression ([56)) does not depend on quantum numbers u> 
and p z . Although result ([55)) has been derived for B' = const field, it can be applicable 
to a spatially slowly varying B' (y) as a good approximation if the gradient variation is 
sufficiently small. 

Let us calculate the total number Af 8 of created pairs with given s defined by (|48|) . To 
do this we go over to the integral 



^ " ■ = L {Lf I dpxdpzdp ° • ' 



Px,Pz,PO 



Taking into account that the exact distribution N n plays the role of a cut-off factor in the 
integral over w, p x , and p z we represent the total number J\f s in the form 

= 2 f dp z N s , p ^ Af s , P z = [ dp x f ff" 2 2 , (57) 

JO (27T) J JO y/UJ^+pj 



where the relation po = 1 + (p z /oj) from ([5]) is used. We obtain the leading contribution 
in ([ST]) as follows 

\/\fj,B'\ ( nm 2 \ / /-— \ 

Af s , Pz = L x L z Tn s ^, n s ^ = ^ exp I - —— j ^V^m^+Pz - \Pz\ ) ■ (58) 

From (|58p . we see that the leading term of density n S:Pz is linear function of the length L y 
for sufficiently small momentum p z , \p z \ <C Umax, that is, the density of the particles created 
per unit space-time volume, n StPz /L y , is uniform. Of course, it is not the case when \p z \ is 
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not small. Thus, we see a complete similarity between the case of a particle creation due 
to a quasiuniform electric field and quasilinear magnetic field for small momenta p z only. 
Using (|57[) . we obtain the total number Af s of created pairs with given s in the form 



Ms 



>/2-l + ]n(l + \/2) 
16?^ 



TL x L z Ly \fiB'f^ 2 exp 



( 



\liB<\ 



2 



) 



(59) 



The total number of created pairs with both s = ±1 is Af =iV +1 + JV_i. 

The vacuum-to-vacuum transition probability defined in (j52")l can be calculated in the 
same way. Then we express it via the total number N as follows 



5 Discussion 

It is quite usual to calculate the probability P v using the Schwinger representation [TO] 



where W is the one-loop effective action. In particular, for the case of of creation of neutral 
fermions with anomalous magnetic moment this approach was used in |12) . We see that both 
the total number N of created pairs and lnP^ 1 , given by (|59")l and (|60p . respectively, are 
finite for the finite space-time volume of field inhomogeneity and independent of magnetic 
field strength Bo. In particular, for uniform magnetic field, B' = 0, both N and lnP^ 1 
are equal to zero. In contrast to the result of Ref. [T2], it follows from our results that 
the arbitrary strong uniform magnetic field is stable with respect of creation of neutral 
fermions with anomalous magnetic moment and this fact does not depend on the space-time 
dimensions. 

We see that in contrast to the case of a constant electric field, the quantities N and 
lnP~ are quadratic in L y . This is a consequence of the fact that the number of states with 
all possible u) and p z excited by the field B' is quadratic in the kinetic momentum \fiB'\ L y . 
That is also the reason why the density of created pairs and the density of Im W per unit of 
length L y are not constant. In this case the divergence of the effective action W as L y — > oo 
is not linear and it is quite difficult to invent in the framework of the Schwinger approach 
a reliable method of regularization of W for a linearly growing magnetic field. We believe 
that this was the main reason of erroneous results in 12J. Note that our approach in the 
framework of QFT can be used to separate the divergent term of Im W as L y — > oo and 
relate it to pair creation, cf., [27] , 

It should be noted that the particle creation in the linearly growing magnetic field rep- 
resents a wide class of physical situations where the magnetic fields have quasilinear hetero- 
geneity in a big enough but restricted areas. One can also see that the leading contribution 
to differential mean number of created pairs in such fields do not depend on asymptotic 
behavior of the magnetic field as the size of the heterogeneity tends to the infinity. 

Particles with opposite values of spin polarization s have opposite accelerations, the same 
is valid for antiparticles. The cases s = +1 and s = — 1 differs only by opposite directions of 
all the motions, and respectively by the opposite dispositions of all the asymptotic ranges. 
Then, unlike the case of particle creation due to the electric potential step, the neutral 
particle (antiparticle) created with different s form fluxes aimed in opposite directions. 
Probabilities of all the processes are equal in both the cases. We see that flux of neutral 
pair created is formed from fluxes of particle and antiparticle of the same spin polarization 
and equal intensity. It is the typical property of neutral particle creation by inhomogeneous 
magnetic field that can be used to observe the effect in astrophysics. 




(60) 



Pv = 
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